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Abstract 

The Cosmic Defect (CD) theory is reviewed and used to fit the data 
for the accelerated expansion of the universe, obtained from the apparent 
luminosity of 192 Snia's. The fit from CD is compared with the one 
obtained by means of ACDM. The results from both theories are in good 
agreement and the fits are satisfactory. The correspondence between both 
approaches is discussed and interpreted. 

1 Introduction 



As it is well known, an extremely important finding of the last decade has 
been the accelerated expansion of the universe. This was rather a surprise, 
mainly based on the observation of luminosity distance of type la supernovae 
(Snia) [1]. Nowadays, the picture which seems to emerge from the data is 
that of an universe which has undergone a transition from a decelerated to an 
accelerated phase, with a relatively recent turning point located at ztr — 0.46 [2]. 
This framework seems to be confirmed by cross-comparison with other pieces of 
evidence [3] . The discovery gave start to an active search for an explanation on 
the theoretical side, within and outside general relativity (GR). An immediate 
effect was to revive the old cosmological constant, A [4]; afterwards, a number 
of evolutionary sons of A or new exotic fields were elaborated, mostly based on 
the idea of "dark energy" [5] [6] [7] [8] [9] . Also various possibilities of alternative, 
modified or extended versions of GR have actively been explored [10]. 

Here, our purpose is to review the existing observational data and some pro- 
posed fits, comparing them one with another and with the results of a recently 
introduced four- vector theory, that we shall call the " cosmic defect" theory, CD 
for short [11]. The CD theory, which has also correspondences in the group of 
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the so called vector "ether" theories [12], will also be revised and recast in the 
following. 

Whenever a theory is contrasted with the data from experience (here, from 

observation) one has to face a mimber of different problems. First of all. comes 
the reliability and cleanness of the data: we shall not elaborate on this, assum- 
ing the discussion to have been effectively conducted in the literature [13]. A 
second subtle issiie is that, even in presenting apparently raw data, underlying 
assumptions often exist, originating in one or another theoretical view: as far as 
possible, we shall try to express the existing information in a model-independent 
way. Finally any theory usually has (a mimber of) free parameters to adjust, in 
order to fit the experiment; of course the more parameters you have, the more 
you will be able to reproduce a given empirical trend, however any choice must 
be cheeked for consistency in as many different physical situations as possible. 

As we shall see, the CD theory gives a reasonably good fit for the Snia 
data, making use of a limited number of parameters and, in the same time, 
offers an interpretation paradigm based on correspondences with known physical 
phenomena without the need of calling in new dark entities. 



2 Luminosity distance, magnitude and redshift 

In the framework of the supernovae observations, a key role is played by the 
concept of luminosity distance, di , which is defined as 



Lobs is the absolute luminosity of the source (released energy per unit time) 
corresponding to the z value measured by the observer, $ is the energy flux 
density (energy per unit time and surface) measured at the observer's site. In 
an expanding universe both energy and time are affected by the expansion so 
that the effective luminosity for the observer, in terms of the absolute luminosity 
at the source is [14]: 

T - 

{l + zf 

In an universe endowed with the typical Robertson Walker (RW) symmetries 
(1) becomes 

di = aors{l + z). 

where oq is the scale parameter at the observer, and rs is the coordinate distance 
of the source from the observer. The latter, written in terms of the distance 
travelled by a light ray, is in turn 



rs = c 



ts W)' 



2 



where of course t is the cosmic time. In terms of the redshift and the scale factor 
we may also write 

da cdz 

where the dot denotes the derivative with respect to t and H = a /a is the 
Hubble parameter. 

It is then easily seen that the luminosity distance is 



di=c{i+z) f (i+Q^=c{i+z) r 

Jo <HU Jo 
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Usually astronomical objects are classified in terms of their magnitude m, 

rather than their luminosity. By definition, the bolomctric magnitude (inte- 
grated over all frequencies) depends logarithmically on the luminosity distance, 
according to the formula: 



m-Ms = 25 + 51ogdi = 25 + 51og |^aoc(l + 2;)y — ) 
25 + 51og {c{l + z) 



'ts 

dC 
HiC) 



(3) 



25 + 5 log c 



(1+^) r dc \ 
m) 



where distances are expressed in Mpc and it is Ho = H (0) and E {z) = 
H {z) /Hq; m — Ms is usually called the "distance modulus". 

The integral in (3) depends of course on the model one uses to describe the 
cosmic expansion. For a dust filled universe in a typical Priedman-Robertson- 
Walker (FRW) scenario it is indeed 



a{t) = anV^T^G^t''/\ (4) 

being p^Q the present matter energy density and G the gravitation constant. 
As a consequence one expects 



(m - Ms)jrjiw = 25 + 5 log 



3c 



{1+z- ym) 



(5) 



If one considers a A-cold-dark-matter universe {AC DM), i.e. an FRW uni- 
verse with a cosmological constant A, it is 

E (z) = ^^l^{l + zf + (6) 

where = Pm/pc represents the ratio between the matter density and the 
critical density (ensuring the flatness of space). The difference Qa = 1 — 
allows for the effect of the cosmological constant. 
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Formula (6) is a special case of the more general 



Eiz) = + zf'^-^\ 

allowing for any number of components of the content of the universe, with 
different equations of state. 

3 The Cosmic Defect theory 



The CD theory is based on the presence of a cosmic (four)-vector field in the 
universe. This vector field is interpreted as the strain flux density in a continuum 
with a pointlike defect^ [11]. The pointlike nature of the defect induces the RW 
symmetry; this very symmetry, together with the defect paradigm, implies the 
vector to be "radial", i.e. everywhere parallel to the cosmic time axis, and 
divergence- free, except at the defect. The norm of the vector, %) will coincide 
with the absolute value of its time component and, according to the divergence- 
free feature, will be: 

X = $. (T) 

where Q is a constant and a is the scale factor of the RW metric. 

The other relevant feature of the CD theory is in the choice of the Lagrangian 
for the spacetime containing the defect. This choice is inspired by the corre- 
spondence between the (bidimensional) phase-space of a RW universe and the 
one of a point particle moving through a viscous fluid [11]. Including the pres- 
ence of "matter" (i.e. whatever is not accounted for by spacetime), the action 
integral is 

S = J (Ke-^'^"^''^'' R + Cmatter) V^d^X, (8) 

with K = d^/WwG and d^x = dtdrdOdip. Explicitly introducing the RW sym- 
metry and considering matter in terms of scalar functions, the Lagrangian read 
out of (8) is: 



Que ^ (a^a -|- ad^) -|- Kofa^a^ + who? , (9) 



where Vk is the part of the Lagrangian which is not affected by any variation 
with respect to the metric, and, in the flat fc = case (polar coordinates), equals 
sinO. The presence of matter is represented by two scalar functions, / and h, 
coupling to spacetime through the constants kq and w. The / function accounts 
for a possible coupling to the rate of expansion of the universe, d. 

The second derivative of a with respect to t, appearing in (9), is easily 
eliminated, once the action is integrated by parts, thus giving a final effective 
Lagrangian of the universe 



^Actually the defect could correspond to any singular hypersurface. 
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— r- + o I + Kofa^a? + who: 



(10) 



Prom (10) the Hamiltonian function is readily obtained, 



n = a—-L 
oa 



Kofa — Qne ^ [ + a 



^^+a^ b?-whaf^. (11) 



As usual, H can be interpreted as the energy content of the system described 
by the effective Lagrangian (10), so that in our case it represents the energy 
content of the universe. The Hamiltonian of an isolated system is a conserved 
quantity, since it is identically 



^-a— o— — -— d-— d=0 
dt da dt da da dd 

Prom now on use will be made of a = a/Q so we write 

6 



Kofa — Qne 



or 



^ - who?= W = const. 



(12) 



(13) 



Prom (13) one directly gets the expansion rate equation: 

W + who? 



Actually, if wc want to recover the usual meaning of the matter term in a 
comoving reference frame, we must choose 

Ko = 0, 

so that the expansion rate can be rewritten as 

W + mha^ 



(14) 



In the absence of a defect we should recover the classical FRW model; for this 
reason, it should be 

' ' (15) 
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where pc^ is the energy density of matter. However, under the same condition 
(<3 = X = 0) eq. (14) gives 



Consistency between (15) and (16) then requires 

w = —1/6, 
h = pc^. 



(16) 
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The final formula for the expansion rate of the universe is 

6/ce-V"«(-^+a)- 

Let us now suppose that the cosmic fluid is made of a number of different non- 
interacting components, each with its equation of state in the form 

Pi = WipiC^, 

where Wi are real positive numbers {wi > 0), and Pi is the partial pressure of 
the ith component. 

The conservation laws imply that 



a, 



S{l+Wi) 

Introducing this relation into (17) we have: 



3(1 + ™,) 



The corresponding Hubble parameter is: 



3(l + »i) 

H = - = - = -^ ' 



(1 + ^) 



3/2 



EPio (l + ^)=^"'^-ra 



l + 6(l + z)Vao 



with ro =-r^. 

In the case of dust (w = 0) and radiation (w = 1/3) it is 



\ 



l + so{l + z)-b 



e-(i+^)Vag 



l + 6il + zf/al 



(19) 



(20) 



The adimensional quantity Eq = Pro/PmO is the present ratio between the radi- 
ation and the matter energy density in the universe, b is to/pmO- 



4 Observations vs theory 

In order to compare theory and observation we make reference to the same set of 
data used recently by Davis et al. [15] and incorporating supernovae analyzed in 
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Figure 1: Fit of distance modulus observations using a standard dust filled 
Friedman-Robertson- Walker universe. The data are from 192 Snia's as ex- 
plained in the text. Vertical bars represent the experimental uncertainties ( 2 
a). The uncertainty on the redshift parameter z would be unperceivable at the 
scale of the graph. The reduced is 2.1276. 



four different groups: 60 from the ESSENCE (Equation of State: SupErNovae 
trace Cosmic Expansion) project [16], 57 from SNLS (SuperNova Legacy Survey) 
[17], 45 nearby supernovae, 30 detected by the Hubble Space Telescope and 
qualified as "golden" supernovae by Riess et al. [18]. As mentioned in the 
introduction, we shall not enter into the discussion of the elaboration of the 
data, but assume them exactly the way they are published or anyway accessible 
considering them as the best available at the moment. 

Altogether we use the luminosity data from 192 Snia [19] which we try to 
fit with theoretical models. 

First we use (5) and obtain the result shown on fig. 1. 

The direct inspection of the graph shows that the data correspond to sys- 
tematically lower luminosities than the ones given by the FRW model, whence 
the accelerated expansion interpretation comes. 

The next step will be to test on the data the KCDM model in its simplest 
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version. To that purpose we use (3) and (6). In practice 



m-Ms = At + 51og(l + z) + 51og / ' (21) 

The result, as it is well known, is better than before, since the reduced 
now is = 1.0295 with a best fitting /i = 43.30 ± 0.03, which correspons to 
iJo = 65.6 ± 0.9 km/sxMpc, and = 0.27 ± 0.03, i.e. 27% of ordinary and 
dark matter plus 73% of dark energy (cosmological constant) in a spatially flat 
universe. For the optimization as well as for the determination of the uncertainty 
of the values of the parameters, use has bcxm made, as for the previous FRW 
case, of a multidimensional nonlinear minimization by moans of the MINUIT 
engine [20]^. The optimization is made minimizing the reduced of the fit. 

Finally we test the CD theory. Use is made of (3) and (20) considering dust 
and radiation, so that 



m-Ms = /i+5 log (l + ^;)+5 log ^ ^ 



i + 6(i + c)Vc^ES 



(l + C)'(l + eo(l + 0-&) 



■dC (22) 



We could treat /x, ao and 6 as optimization parameters, however the value 
to be introduced for eo is the one currently agreed upon, excluding any dark 
contribution: Eq ^ 10^^. Of course as far as z is in the order of a few units (as 
it is the case for Snia's) the radiation term in the denominator of the integrand 
is negligible, so that also the contribution of h may be embedded in \x and the 
free parameters remain /x and only. The result of the optimization process 
is = 43.26 ± 0.03 and ao = 1.79 ± 0.04; the reduced is = 1.092, almost 
as good as for KCDM. The graph is shown in fig. 2. 



4.1 Effective gravitational coupling constant 

The meaning of the choice of Snia's for cosmological analyses relies on the 
assumption that they are "standard candles", which means that the relation 
between their light curve and luminosity is thought to be always the same, now 
as well as at epochs corresponding to high cosmic redshifts (actually, at the 
moment, the highest redshift for an Snia corresponds to 2; ~ 1.8). If it is so, in 
an expanding universe we expect any standard candle (as an SnIa is assumed 
to be) to appear dimmer and redder according to the peculiar expansion law. 

Indeed the currently assumed explosion mechanism for la-type supernovae 
makes their peak luminosity be proportional to the mass of nickel synthesized, 
which, to a good approximation, is a fixed fraction of the Chandrasekhar mass 
[21]. The Chandrasekhar mass, in its turn, is proportional to G^'^l^ [22]. In 

^The open source routine we used, due to G. AUodi of the university of Parma, 
is named fmiiniit, is called from within MATLAB, and may be retrieved from 
ftp://ft p. fis.unipr. it/pub/mat lab/fminuit. mex 
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Figure 2: Fit of distance modulus observations using the CD theory with two 
free parameters. Radiation is overlooked with respect to dust. Symbols as in 
Fig. 1. The reduced is 1-092. 



an expanding universe it is reasonable to expect that an effective value of G, 
somehow accounting for the expansion, must replace the usual constant in the 
formulae we use to describe locally the behaviour of matter under gravity. If 
it is so a consequence is that the absolute magnitude of a source at a given z 
differs from the one for z — Q. In practice, including the fact that the extension 
of the light curve (characteristic time) of an Snia is approximately proportional 
to the square root of the Chandrasekhar limit [23] , it is 

M (0)^ - Ms (^) = ^ log 
and the distance modulus becomes: 

15 G iz) 
m{z)-Ms = 25 + 5 log + — log— 

4 Go 

For the above reason it is important, in any model we use in order to describe 
the evolution of the universe, to determine the explicit form of G{z). This 
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is indeed a delicate task, since everything has to be worked out in terms of 
appropriate approximations. Attempts have been made to solve this problem 
in general terms for any scalar-tensor theory leading to acceptable results for 
limited ranges of z values [24]. We would like to follow a different path. 

Actually at the scale of a stellar system the effect of the cosmic curvature 
appears only as a very tiny perturbation of the usual stationary state solutions 
of general relativity. Considering for instance the example of a spherically dis- 
tributed bunch of matter in a background RW universe, we expect that locally 
the induced metric be essentially 

ds^ = (1 - / (r, I)) di^ - a" (t) [^ ^X{ri) + ''^'^^^ + ^^"^ ^^'^') ^^^^ 

with an extremely weak dependence of / and h on time also. If tq is the given 
cosmic time we use t = t — tq as the time variable now. Far away from the 
local source it must be /, /i ^ so that the pure RW metric is recovered. On 
the other side if it were a — ao ~ constant we could incorporate its value in the 
definition of r and from the field equations we would obtain the Schwarzschild 
solution f — h — 2Gm/c^r. 

As written above, ao may be absorbed into a rescaling of r and similarly 
Go may be included into / and h. Now we may think to develop a, /, and h 
in a power series of t; if we do so, including the condition /, h << 1, we end 
up [25] with the interesting result that the zero order (in the typical evolution 
times of the gravitational phenomena we wish to describe) is independent from 
such parameters as Hubble's. In practice rapid events, such as the accretion of 
matter preluding to a gravitational collapse, always happen to be controlled by 
the universal value of G and no correction is expected on the distance modulus. 
Any effect of the expansion of the universe (differentiating one model from 
another) emerges only for times comparable to the cosmic evolution times. 



5 The Hubble parameter and the age of the uni- 
verse 

Reconsidering now the explicit spelling of the parameters appearing in the CD 
theory used to draw fig. 2 we see that it is 

^ = 25 - 5 log c + I log ■ 



2 ''Pmoil-b) 
whence 

Pmo{l-b) = fl0-t(''-2^) (24) 



= (8. 5 T 0.2) X 10-27 kg/in 

The "visible" matter density in the universe is commonly assumed to be around 
~ 10-2'^ 10"^^ kg/m^ which means that b must be -10. 
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Then, introducing (24) into (20) and evaluating for ^ = we obtain 



^ ^ I PmO (1 + gp - b) 
° ~ Y 6Ke-i/«o(l + 6/a6) 
= (62. 8 T 1-7) km/(s x Mpc) 

which is an acceptable resuh (eo has been neglected with respect to b). The 
corresponding Hubble time is 15.6 Gy. 

Of course we should determine the age of the universe using the CD model; 
this can be done by means of (18) by integration: 



= (9.0 ±0.2) Gy 

The final numerical result has been obtained neglecting eoao with respect to 
the other terms in the denominator of the integrand. The value falls rather 

short as compared to the age of globular clusters, which fact may probably be 
interpreted as an inadequacy of the model at very early cosmic times. 



6 Conclusion and discussion 

We have fitted the apparent luminosity data from Snla's with the values pre- 
dicted by the ACDM and the CD theories, comparing both with a traditional 
FRW universe. The result is of course partly known, but we see now that also 
CD improves with respect to FRW and gives a fit comparable with the one of 
ACDM. Using the same data and the same number of parameters we obtained 
similar values of the riduced x^'s suggesting the idea that CD also is a viable 
theory. It is however true that the apparently small difference of the reduced 
X^'s of the fits corresponds to a rather big difference in the full that, when 
analyzed on the light of statistical information criteria, such as the Akaike Infor- 
mation Criterion (AlC) [26] and the Bayesian Information Criterion (BIC) [27], 
enhances the distance between the two theories in favor of ACDM. In the same 
time it is also true that both reduced values of are bigger than 1; furthermore 
the Hq values obtained from observation using different methods are systemati- 
cally higher than the ones of the 2-parameters best fits above. The most recent 
data from WMAP [28] yield Hq = 73.2tf^ km/sxMpc which is consistent with 
a number of other results produced by different methods and indicators (like 
Snl, Snll, Cepheids in nearby galaxies, Sunyaev-Zeldovitch effect, X rays from 
clusters, gravitationally lensed systems) all quoted in [28]. The central values 
from these different observations range from 72 to 76 km/sxMpc and in general 
the historical evolution of the estimated values of the Hubble constant seems 
to progressively converge towards something around 75 km/sxMpc^, which is 

^Look for instance at http://cfa-www.harvaxd.edu/''huchra/hubble/ 
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^ 15% more than the resuhs got by means of the fits in this paper. If the "ex- 
perimental" value of Hq were used in the fits (so reduced to 1-parameter ones) 
the agreement with the data would consistently worsen both for ACDM and for 
CD. In practice there is something missing beyond the details of the theories 
and their interpretation, which deserves investigation and insight. 

The ACDM is indeed different from the CD theory: the former assumes in 
the universe the presence of a cosmological constant corresponding to a sort 
of uniformly and homogeneously distributed dark energy; the latter interprets 
space-time as a continuum with a cosmic defect inducing a strained state con- 
taining both the symmetry and the non-uniform expansion rate. Besides this, 
we know that ACDM requires also that the matter content in the universe be 
one order of magnitude bigger than what expected from baryonic particles only. 
In the case of the CD theory, instead, we saw that the ordinary matter density 
is combined with the effect induced by the defect via the b parameter (see 24) , 
so that, in a sense, it gives rise to an effective matter/energy density one order 
of magnitude bigger than the actual one. Adding the fact that one can inter- 
pret the strained state induced by the cosmic defect as being the equivalent of a 
non-uniform (in time) dark energy, we see that in facts the principle difference 
between ACDM and CD could not be that deep. The CD theory already proved 
to correspond to vector theories developed with different motivations and within 
a different scenario [12]. The real difference between the theories, not consider- 
ing the details, is in the end in the fact that CD tries to give an "explanation" 
for the dark energy, that seems to permeate the universe, within a consistent 
paradigm, which is the one of the defected continuum with its properties. 

Of course this test of the CD theory is limited to the Snla data and conse- 
quently to a limited range of z values, and the poor result obtained for the age 
of the universe seems to indicate an inadequacy of the theory at high redshift 
values, where probably a better treatment of the matter content is in order. 
The result with the type la supernovae is however encouraging. 
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